Local perturbation in a Tomonaga-Luttinger liquid at g = 1/2: 
orthogonality catastrophe, Fermi-edge singularity, and local density of states 
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The orthogonality catastrophe in a Tomonaga-Luttinger liquid with an impurity is reexamined 
for the case when the interaction parameter or the dimensionless conductance is g = 1/2. By 
transforming bosons back to fermions, the Hamiltonian is reduced to a quadratic form, which allows 
for explicit calculation of the overlap integral and the local density of states at the defect site. The 
exponent of the orthogonality catastrophe due to a backward scattering center is found to be 1/8, 
in agreement with previous studies using different approaches. The time-dependence of the core- 
hole Green's function is computed numerically, which shows a clear crossover from a non-universal 
short-time behavior to a universal long-time behavior. The local density of states vanishes linearly 
in the low-energy limit aX g = 1/2. 
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I. INTRODUCTION 

One-dimensional interacting j-, fermion systems, 
Tomonaga-Luttinger (TL) liquidsJTH have recently at- 
tracted much attention due to their anomalous resnonse 
to local perturbations. Recent extensive studiescFLl on 
transport properties of TL liquids with an impurity re- 
vealed that repulsively interacting fermions have vanish- 
ing transmission probability through a potential barrier 
in the low-energy limit. This is because the interac- 
tion between fermions strongly enhances the backward 
scattering at the barrier. Thus, a single defect effec- 
tively cuts a TL -liquid into two disconnected ones at 
zero temperaturcLl This implies that the local density of 
states (LDOS) at the defect is reduced for low energy, 
and according to Kane and FisherQ it shows a power-law 
energy dependence. 
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where g is a parameter characterizing the TL liquid. This 
picture was^-however, questioned recently by Oreg and 
Finkel'stein,E3 who claimed based on a mapping to a 
Coulomb gas problem that the LDOS at the defect is en- 
hanced, rather than suppressed, in the low-energy limit 
for weakly interacting fermions. This controversy moti- 
vates us to reexamine this issue. |— . 

The orthogonality catastropheO in a TL liquid is an- 
other interestiag-subject which has been discussed by sev- 
eral authors .tHrii3 They showed that the overlap between 
the ground state of a pure TL liquid |p) and that of a TL 
liquid with a single scatterer |s) vanishes in the limit of 
large system size: 
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where L is the length of the system. The exponent ■^p 
is due to the f^mwd-scattering potential and depends 
on its strength.EJO It can be calculated directly using 



a unitary transformation. The other exponent jb due 
to the backward scattering is believed to be indepen- 
dent of the stpength of the potential and take a univer- 
sal value, l/8lli"ll3 In Refs. |l|,|l|, and |l| the exponent 
7b is calculated by assuming that a backward scatter- 
ing center can be replaced with aiuimpenetrable poten- 
tial barrier. Oreg and Finkcrstein,E3 however, questioned 
the validity of the assumption and argued that the ex- 
ponent of the Fermi-edge singularity due to a backward 
scattering center is zerch-which implies 7b = 0. On the 
other hand, Kane et a/.ll3 used a renormalization-group 
equation which becomes exact in the limit of weak repul- 
sive interaction between fermions. They could describe a 
crossover from the high-energy regime to the low-energy 
regime, and obtained the same exponent 7b — 1/8 in the 
low-energy limit. The result of a repent direct numeri- 
cal calculation of the overlap integraO is also consistent 
with 7b = 1/8. 

It is known that, when the TL-liquid parameter g is 
1/2, the bosonized Hamiltonian containing a nonlinear 
term representing the backward scattering can be_t;rans- 
formed to a quadratic Hamiltonian of fermions.t3 This 
is essentially the same technique as the Ernery-Kivelson 
solution of the two-channel iCondo problemHj The exact 
results pQir the conductanceu and non-equilibrium noise 
spectraEJ were obtained using this refermionization tech- 
nique. It is thus natural to expect that exact calculation 
should also be possible for the above-mentioned prob- 
lems. The purpose of this paper is to show that this is 
indeed the case. 

The structure of this paper is as follows. After in- 
troducing a model of interacting fermions in Sec. ||, we 
discuss in Sec. HI the exact low-energy behavior of the 
LDOS for g = 1/2. For g ^ 1/2 we show that Eq. (|) fol- 
lows from the assumption that the phase field is pinned 
at the defect site. The importance of zero-modes is em- 
phasized. In Sec. IV we calculate 7b analytically for 



g — 1/2 without assuming the nature of the low-energy 



fixed point. We find jb = 1/8. The so-called core- hole 
Green's function is then computed numerically in Sec. |y|, 
which shows a clear crossover from short-time to long- 
time regimes. We show in Sec. VI that the exponent of 
the Fermi-edge singularity due to backward scatteri ng is 



also given by 75. We summarize the results in Sec. VII 



II. MODEL 

In this section we introduce a model of interacting spin- 
less fermions and briefly explain the bosonization rule to 
fix the notation. 

The Hamiltonian of our model is given by 
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where ipLiB.) describes left-going (right-going) fermions, 
: A : represents normal-ordered operator A, and Xp 
(Ase*^) is the forward-scattering (backward-scattering), 
potential. Following the standard bosonization rule,E^ 
we express fermions ip^, in terms of bosonic operators: 
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where a is a short-distance cutoff. The bosonic fields 
satisfy the commutation relations [ipL{x),ipL{y)] — 
-i7TSgn{x - y), [ipR{x),(pR{y)] = i7rsgn(x - y), and 
[(pL{x),ipR{y)] = 0. The operator ?7^t's are Majorana 
fernfions corresponding to zero modes of bosons, which 
are needed to ensure the anticommutation relation be- 
tween ipL and 'ipB.- They satisfy {rjL, tjr} = and rjf^ = 1. 



We then introduce new bosonic fields as 
1 
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[iPr{x) + ipLix)], 
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which obey [0(x),n(j/)] = i5{x — y). With these fields 
the Hamiltonian can be transformed to a bosonic form. 
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We then introduce another set of bosonic fields ip± (x) 







The parameter g is related to 52 and g^hy g = [(1 + g^ — 
52)7(1 + 54 + 52)]^'^ with (ji = gi/2TTVF- Since the in- 
teraction is repulsive, g is less than 1. The renormalized 
velocity is given by w = vp[{l + g^)"^ — {g2)'^Y^'^- 



f±{x) == -| I f — - V5J [fR{x) ± fL{-x)] + f — + V5 ) [fR{~x) ± ifL{x)] 
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These fields satisfy [<^j^{x), ip-\.{y)] — [ip-{x), ip-{y)\ = — i7rsgn(a; — y) and [(^+(a;), ip-{y)\ = 0. The advantage of using 
Lp± is that we may separate the Hamiltonian into two commuting parts, H = Hp -\- Hb , where 
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The fermion field at a; = may be written as 
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III. LOCAL DENSITY OF STATES AT A SCATTERING CENTER 



In this section we calculate the following correlation function: 



D{t)^{geW"'^\Q)e-^"'i,mm). 
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where \ge) is a ground state of H. The LDOS is given by p{ijj) — j {dto / 2n)e'^'^'^ D{t) . In general we expect D{t) ex 
g-iAt^-i^ for i — > oo. Since H has gapless excitations, we know that A must be zero. Thus, we will not pay attention 
to A and concentrate only on the exponent v in the following discussion. 

Since Hp and Hb commute, the correlation function is factorized into two parts as D[t) = -:^^D p{t)D sit) , where 
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Here $_ = lp^{0)/^/Ig, <f>+ == y^gj2ip+{0), and |F) (|B)) is a ground states of Hp (Hb). The Hamiltonian Hp is 
related to a free Hamiltonian by a unitary transformation as UHpU' = Hp + const, where 
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and 



U — exp 
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This means |F) = W\Fq) with |Fo) being the ground state of Hp . We thus get 
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As pointed out in Ref. [lO|, t he forward-scattering potential does not affect the LDOS. 
Next we rewrite Eq. ([12b|) as 
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where Hb = VlHb^l — Hb{)^b — * —As). Note that this 
sign change of the cosine term is a direct consequence 
of the anticommutation relation {V'LiV'-r} — 0. At this 
point we may set rj^ripi = —i because only the terms 
involving even powers of rjLVR will contribute to DB(t) 
when Eq. ( p^ ) is calculated perturbatively in powers of 

\b- " 

obtain 



Ar. We then shift (p+{x) -^ <(5+(a:) + :^(f — 0) and 
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where 
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and we have used the fact that the ground state of iJ+ , 
1+), is invariant under if^ — + —ip+. It is useful to trans- 
form Eq. (|l7|) further to the form 
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We first consider the case of 5 = 1/2. A cru- 
cial point in this case is that the cosine term becomes 
e'v+io) _^g-*v+(o) Therefore, fermionizing the chiral bo- 
son Lp^ as 



o«V+(a:) 
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where 
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where ry is a Majorana fermion, satisfying 77^ = 1. This 
leads to a simple relation, r/H+g = i7_. Note that 
iJ+ is a quadjTatic Hamiltonian, which can be easily 
diagonahzed:E3 
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where ^k = vk and ip+ix) = J{dk/^/2n)e ^'^^ak- For 
later convenience we write the transformation rule hereicHI 



dk^k ( cick + didk ) + const, 
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where A: > 0, F = A^/(7raf), and Ck and dk satisfy the ordinary anticommutation relation. The ground state 
the vacuum of Ck and dk ■ 

Using Eq. (vW, we rewrite Eq. (M) in a fcrmionic form, 
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From Eqs. (24a) and (24b), the second term becomes 



-const. 
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which is an irrelevant operator with scaling dimension 
2. To find the long-time behavior of DB{t), we can thus 
treat Eq. (|2^) as a small perturbation. The lowest-order 
calculation then gives, for Tt ^ 1, 
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For g ^ 1/2 (0 < g < 1) we take a different approach. 
We assume from the outset that the phase field (p+ is 
pinned at a: = by the cosine potential in H+ (18), as 
in Rcfs. |lj,|l^, and |l^. Wc thus replace the cosine by a 
term which is easier to deal with. A convenient choice is 



(28) 



Hm — 



Note that the 1/t-dependence of the first term comes 
from the correlator (+|?7(0'7(0)|+), which-i also ap- 
peared in the two-channel Kondo problcmB3 Combin- 
ing Eqs. (|15|) and (|2|), we get D{t) = -2/{Tr^vTt^) for 
Ft 3> 1, which implies 
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where M should be a characteristic energy scale at which 
the cosine term becomes of the order of the band width 
(M = F for g = 1/2). It immediately follows from the 
scaling equation dXs/dl = (1 — 5)A_b that 
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for cj <C F. This is consistent with Eq. (|lj). We see 
that the single scatterer at a; = indeed depletes the 
low-energy excitations around it. 



Since H^ is a quadratic Hamiltonian, it is easily diago- 
nalized as Hm = J dk(_k V^k^k + Pk^A "^'^^^ 



dk 
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(32) 



and n = {l/gv)dll/dt, where ak and (3k satisfy the ordinary commutation relations of bosons. The phase shift is 
given by Sk = tan~^{gM/2vk). Note that 5k -^ 7r/2 as fc — > 0. 
Let us denote the ground state of Hm by |0m)- We then find 

24 
(0M|5,^+(0,t)5,^+(0,0)|0M) = 27rg(0Af|n(0,i)n(0,0)|0M) = 2M2y2t4 (33) 



for Mt :$> 1, implying that dx^p+iO) is an irrelevant op- 
erator with dimension 2. This is consistent with the ob- 
servation made in Eq. ([27|). In fact, this is an expected 
result because ip+ is pinned at a; = 0. We may thus 
use H- instead of H- to obtain the long-time asymp- 
totic behavior of Db^I) in Eq. dl9). It is also important 
to note that e'*+ is not fluctuating too much and can 
be regarded essentially as a constant because (^+(0) is 
pinned. In fact, we find 



(OM|e**+|OM>-(OAf|exp iVW5</.(0) |0m 



e-'gaM 
2v 



(34) 



for aM <C v, where 7 = 0.577... is Euler's con- 
stant. Note that, at g = 1/2, we get (+|e**+|+> == 
— (Ab/ttd) ln(u/ar), which is consistent with Eqs. ( |3l| ) 
and (|34|). Hence, from Eq. (|l9|), we get 
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where T^ is a unitary operator which shifts <j){x) — > 
4>{x) + ^^^. The rhs of Eq. (|35| ) is known to decay as 
~ i~i/2ff£5| This result can be easily obtained using the 
following representation for V: 



V = exp 
= exp 



dxH{x) 
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We therefore conclude D{t) ex t^^/^, from which Eq. (|l|) 
follows. We emphasize that the above calculation should 
give the exact value of the exponent, although the am- 
plitude may not be correct. 



IV. ORTHOGONALITY CATASTROPHE 



In this section we discuss the orthogonality catastro- 
phe for the special case oi g = 1/2. We calculate the 
overlap integral |(p|s)|2 = |(Fo|F)|2 x |(0|+)p, where |0) 
is the ground state of the Hamiltonian Hq = H^Ix^^q. 
It is almost trivial tpJiajd 7f in Eq. (||) because {Fo\F) = 
(Fo|[/t|Fo). WegetBia 



^2af'- 



IF = 2g 



V 27rw 



(37) 



Hence our problem is reduced to calculate the overlap 
(0|+). In the fermion language, Hq is 



Hr, = iv 



dxipl{x)—ilj+{x), 



(38) 



and |0) is the filled Fermi sea. Then the ground state of 
Hj- can be written as 



K) = Texp 



e''*H'{t)dt 



|0), 



(39) 



where e is positive infinitesimal and H'{t) = e*^"*(i7+ — 
i?o)e *^°*- Using the linked-cluster theorem, we can 
write the overlap integral as 

(0|+)=cxp[Gc(0,-cx3)], (40) 

where Gc(0, — cxd) is a sum of connected ring diagrams. 
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Here A = Xb /V^na and the propagators so{t) and go{t) are given by 
so(i) = (0|T,7(t),7(0)|0)-sgn(t), 
go{t) - (0|T[^+(x- 0,t) - V|(0,t)][V'+(0,0) - V^|(0,0)]|0> = 



Trv[t — i£Sgn(t)] ' 
where e is positive infinitesimal. Differentiating Eq. (|4l|) with respect to A, we obtain 
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where g{ti, ^2) is a solution of a Dyson equation 
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Since Eq. (53) contains double integral, working in real time is not so convenient as it is in the Fermi-liquid case.cJ 
On the other hand, the Fourier transform of Eq. (H) contains only a single integral: 
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This equation can be solved in the hmit e ^ in the standard wayl 
We first introduce functions g± by 
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We can then express Eq. (Bsh as 
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A solution of this equation with correct analytic properties is 
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where 5 is positive infinitesimal and 
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With this solution Eq. (ES) becomes 
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where t = i2 — ti and we have integrated over (ti + t'i)l2. As pointed out by Hamann.Ea in the next step in which 
we perforin the w integral, it is important to keep e finite while taking the limit b — > +0: 
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After replacing v/{\v\ ~ iV) by v\{\v\ - iVy^ - {v - t^Y^\ and ln(l - i^\ by ln[(l - i^^l (l + ^7)], we integrate 



over T to obtain 
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where we have introduced the high-energy cutoff A ~ vj a aHjd the , low-energy cutoff E^ ~ vjL. From Eqs. ( [40[) and 
( ^ ) we get 7b = 1/8 in agreement with the previous studies .IlirliZlllj Note that the quantity Eq = — (r/2p4[ln(A/r) + l] 
appearing in the first term is equal to the difference between the ground state energies of iJ+ and iJo -tHl 

Since b{E) = tan^^(r/_E) in Eq. (Ell) is the phase shift for fictitious chiral fermions due to the cpuxiling Xb in 
Eq. (P4), the above calculation impHes that jb = ^['^(O)/''']^! i"^ contrast to the Fermi-Hquid resultOO 7Fermi — 
[(5(0)/7r]2. The extra factor 1/2 in our result can be traced back to the peculiar form of the scattering term in 
Eq. (E2). Only the combination V'+ — V'+ interacts with 77, and the other combination ■0+ + ■i/;'|_ is decoupled. Hence 
only half oi the degrees of freedom have the phase shift (<5(0) = 7r/2), giving the factor 1/2. 

As pointed out by Matveev,E2l the. Hamiltonian (pa) is equivalent to the effective Hamiltonian of the two-channel 
Kondo model in the Toulouse limit J£3 where the Majorana fermion 77 corresponds to the xy-component of the impurity 
spin. Thus our calculation also applies to the orthogonality catastrophe in the two-channel Kondo problem in which 
J± is turned on and off while Jz kept constant. 



V. CORE-HOLE GREEN'S FUNCTION 



Next we calculate the core-hole Green's function, 

G{t) = (0|e*^"*e-'^+*|0) 
for g = 1/2. Using the linked-cluster theorem again, we get G{t) — exp[Gc(t, 0)], where Gc(i,0) is 
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This time we differentiate Eq. (|53| ) with respect to t to 
get 



-G,{t, 0) = A^ y dhgoit - ii)s(ii), (54) 



where s(ii) is defined for < ti < t and is a solution of 
a Dyson equation. 



^(^i) = -i-:^PMi3Mi. 
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From this equation we can easily show that s{ti) 
s{t - ti) and s(+0) = -1. Thus Eq. (Hj) becomes 



io.m^'.-L ['"<.■'"> 



dt 



2m 



h 



(56) 



Here the first term comes from the real part of go in 
Eq. ( pb| ). 

For short times Ft ^ 1, we can solve Eq. ( pq ) pertur- 
batively. Up to order (Ft)'^ we obtain 



Gc{t,0) = i 



rt 

2^ 



In 



iFt+l(Ft)2, (57) 



where tc is a short-time cutoff ^ 1/A. This expansion, 
however, starts to fail around Fi ~ 1. From the analysis 
in Sec. |rv|, for F,tT::S2L-l we expect Gc{t,0) to approach 



-iEot - I In(Fi) 

The crossover from the short-time to the long-time 
regimes can be seen most conveniently by solving Eq. (|55| ) 
numerically and putting the solution into Eq. (p9). Note 
that the integral in Eq. (jSq ) is well-defined because 
Inis(ti) ^ ii|lnti| for ti -^ 0. Figure 1 shows the t- 
dependence of the real part of {d/dt)Gc{t,0) computed 
in this way. ft clearly exhibits the crossover at Fi ^ 1 
from the short-time behavior, Eq. (pq), to the long-time 
asymptote, Re[dGc{t,0)/dt] = -l/8t. 
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FIG. 1. Time evolution of the core-hole Green's function. 
There is a clear crossover at Pi ~ 1. The dashed line repre- 
sents Re[dGc/drt] = -l/(8ri). 



VI. FERMI-EDGE SINGULARITY 

In this section we briefly discuss the Fermi-edge singu- 
larity for g < 1 to show that the exponents can be easily 
obtained from the analysis of Sees. IV and ^. Here we 
are concerned with the correlation function 



I{t) = (.9o|e^(^- +^")V(0)e-^(^-+^-)*7A^(0)|.9o), (58) 
where 



Sec.|^ 

where 



Si 



go) = l-fo) <Xi |0). Following the same path as in 
Bie write the correlator as I{t) ~ ^^lF{t)lB{t), 



Ipit) = {Fo\e'"- 



.iW"'t„-i<i?- 



Ue 



-iH'-"'t 



uh 



i*_ 



Fn 



(59) 



^i*l^^^=(^ + ^^^)'and 



/s(t)-2(0| 



,iHot„-iH-t 



|0) 



-2cos6'(0|e*-"°*e-'^-*e'^**+|0). 



(60) 



We expect that /^(i) should decay as Isit) (x t "^^ in 
the long-time limit. We now notice that the first term 



in Eq. ( pQ ) is similar to the core-hole Green's function 
discussed in Sec. 0. As we saw in Fig. 1, it should decay 
as ~ t~^ with 7 being the exponent of the orthogonal- 
ity catastrophe between |0) and the ground state of i/_: 
|(0|— )p oc L~^ . The latter state has a finite overlap with 
the ground state of i/_, because dx^fiiO) [oc (iJ_ — H^)] 
is an irrelevant operator around the fixed point of H^. 
This means 7 = 7s = 1/8. Since the second term in 
Eq. ( pO| ) contains extra factor, e^**+, at least it is not 
larger than the first term. Hence we conclude vb = 1/8, 
in agreement with Refs. O and ll^. The fact that vb 
equals 7^ is a direct consequence of the pinning of ip^ at 
a; = 0. Therefore the insertion of ip^ part of the fcrmion 
field, e**+, does not change the exponent. On the other 
hand, vp is not equal to "fp because the forward scatter- 
ing potential is a marginal operator. 



VII. CONCLUSION 



In this paper we have studied the low-energy behav- 
ior of the LDOS at the location of a scattering center 
and the orthogonality catastrophe due to a sudden local 
perturbation. The characteristic, anomalous low-energy 
(long-time) properties were obtained by exact calcula- 
tions for g = 1/2 by mapping the bosonized Hamiltonian 
back to a fermionic quadratic Hamiltonian. This method 
has allowed us to describe the crossover from the weak- 
coupling (short-time) to the strong-coupling (long-time) 
regimes. The exact results obtained for g = 1/2 agree 
with the previous studies based on the assumption that 
the phase fields are completely pinned at the impurity 
site in the low-energy limit. The agreement implies that, 
to describe the low-energy physics, it is sufficient to use 
an effective model which incorporates the perfect reflec- 
tion by the local potential. We conclude that jb — 1/8 
and p{llj) (x lu'H^ for g < 1. It seems that the mapping 
to a Coulomb gas problem used in Refs. EO and Iq makes 
it difficult to capture the Majorana fermions which have 
played an essential role in this paper. 

After completion of thiSjSSiork the author became aware 
that Fabrizio and Gogolinc3 obtained a similar result on 
the low-energy behavior of the LDOS, Eq. (Q). 
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